Statics and Dynamics of Yukawa Cluster Crystals on Ordered Substrates 
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We examine the statics and dynamics of particles with repulsive Yukawa interactions in the 
presence of a two-dimensional triangular substrate for fillings of up to twelve particles per potential 
minimum. We term the ordered states Yukawa cluster crystals and show that they are distinct 
from the colloidal molecular crystal states found at low fillings. As a function of substrate and 
interaction strength at fixed particle density we find a series of novel crystalline states that we 
characterize using the structure factor. For fillings greater than four, shell and ring structures form 
at each potential minimum and can exhibit sample-wide orientational order. A disordered state can 
appear between ordered states as the substrate strength varies. Under an external drive, the onsets 
of different orderings produce clear changes in the critical depinning force, including a peak effect 
phenomenon that has generally only previously been observed in systems with random substrates. 
We also find a rich variety of dynamic ordering transitions that can be observed via changes in 
the structure factor and features in the velocity-force curves. The dynamical states encompass 
a variety of moving structures including one-dimensional stripes, smectic ordering, polycrystalline 
states, triangular lattices, and symmetry locking states. Despite the complexity of the system, we 
identify several generic features of the dynamical phase transitions which we map out in a series of 
phase diagrams. Our results have implications for the structure and depinning of colloids on periodic 
substrates, vortices in superconductors and Bose-Einstein condensates, Wigner crystals, and dusty 
plasmas. 

PACS numbers: 82.70.Dd,64.70.Rh,64.60.Cn,64.70.Nd 



I. INTRODUCTION 

The creation of new types of crystalline or partially 
ordered states and the dynamics of assemblies of inter- 
acting particles have attracted much attention both in 
terms of the basic science of self-assembly and dynamic 
pattern formation as well as for applications utilizing self- 
assembly processes. One of the most extensively studied 
systems exhibiting this behavior is assemblies of colloidal 
particles, where the equilibrium structures can be tuned 
by changing the directionality of the colloid- colloid in- 
teractions P, Q ■ Since it can be difficult to control and 
tune the exact form of the interaction, another approach 
is to use colloids with well defined interactions that are 
placed on some type of ordered substrate. Optical trap- 
ping techniques are one of the most common methods of 
creating periodic substrates for colloids iSj]. Studies of 
colloidal ordering and melting for one dimensional (ID) 
periodic substrate arrays have revealed ordered colloidal 
crystalline structures as well as smectic type structures 
where the colloids are crystalline in one direction and 
liquidlike in the other |4l-[l3|. These experiments show 
that the substrate strength strongly influences the type 
of colloidal structure that forms and that as the substrate 
strength increases, the resulting enhancement of fluctua- 
tions can induce a transition from crystalline to smectic 
order [1, 0, Q. Numerous different colloidal crystalline 
structures can also appear on ID substrate arrays of fixed 
periodicity when the colloid density is varied [3, 0| . 

More recent studies addressed colloidal ordering on 
two-dimensional (2D) periodic substrates [13 17]. In 
these systems, the filling factor / is defined as the num- 
ber of colloids per substrate minimum. For integer fillings 



f — n, the colloids in each minimum can form an effec- 
tive rigid n-mer, such as a dimer or trimer [l6l - l22j . The 
n-mers have an orientational degree of freedom, and de- 
pending on the effective interaction between neighboring 
n-mers, all the n-mers may align into a ferromagnetically 
ordered state, sit perpendicularly to their neighbors in an 
antiferromagnetically ordered state, or form other orien- 
tationally ordered states. For square and triangular sub- 
strate arrays, n-mer states have been studied up to n = 4 
[l9l42lj : however, it is not known what structures 
would form at higher fillings when the simple picture of 
rigid n-mers no longer applies. Studies of the ordering of 
bidisperse colloidal assemblies with two different charges 
on 2D periodic substrates produced novel ordered phases, 
while a pattern switching could be induced by application 
of an external field [23 |. Similar pattern switching also 
occurs for colloids with monodisperse charges under ex- 
ternal driving [23|. The colloidal n-mer states have been 
termed "colloidal molecular crystals," and for conditions 
under which they loose their orientational ordering, they 
are referred to as "colloidal plastic crystals." Colloidal 
molecular crystals appear for integer fillings f — n. At 
fractional fillings such as / — 3/2, 5/2, or 7/2, it is pos- 
sible for ordered composite states to form containing two 
coexisting species of n-mers; however, for other fractional 
fillings, the srotem is frustrated and the n-mer states are 
disordered [2J]. Other studies have shown that novel 
orderings appear when the 2D substrate array has qua- 
sicrystalline order [2^. 

Once colloidal crystals have formed on a substrate, the 
driven dynamics can be explored by applying an external 
driving field to the sample. A variety of dynamical lock- 
ing phases can occur in which the colloids preferentially 
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flow along symmetry directions of the substrate '26'] . As 
the filling fraction / is varied, a series of peaks in the crit- 
ical force needed to depin the colloids occurs at integer 
values of / indicating the existence of commensurability 
effects [27]. Recent experiments with strongly interact- 
ing colloids on two periodic arrays show that kink-type 
dynamics can occur near / = 1 [281 ] . It would be interest- 
ing to explore higher fillings where new types of dynamics 
could emerge [29| . 

Many of the same types of phenomena found in col- 
loidal molecular crystals can also be realized for other 
systems that can be modeled as interacting particles in 
the presence of a 2D periodic substrate. For example, 
the antiferromagnetic ordering of dimer colloidal molec- 
ular crystals on a square substrate was reproduced using 
vortices in Bose-Einstein condensates confined by opti- 
cal traps with two vortices per trap [s^j as well as with 
vortices trapped by large pinning sites in type-II super- 
conductors ^3 1,-3 3,]. Experimental and numerical studies 
of molecular ordering on periodic substrates show similar 
orderings [s^ ]. Other systems where similar states could 
be realized include classical electrons or dusty plasmas 
with some form of substrate as well as crystalline cold 
atoms on optical lattices. 

To our knowledge, previous studies of Yukawa inter- 
acting colloidal molecular crystals have focused only on 
systems with up to four colloids per trap. For the case 
of three colloids per trap, only a limited number of stud- 
ies have considered the dynamics, and even in this limit 
there are several new features that we describe for the 
first time in this work. We show that at high fillings, 
the rigid n-mer picture breaks down and new cluster and 
ring states form. Several general features of the statics 
and dynamics emerge at these larger fillings which are 
independent of the specific filling. 

One of the key findings in our work is the develop- 
ment of orientationally ordered shell structures at higher 
fillings for the 2D arrays. The development of particle 
shell structures was studied previously for repulsive par- 
ticles in isolated individual traps, in systems that include 
classical charges fssi] , Wigner islands [11] , dusty plasmas 
[37j , colloids ^38r;42j , charged balls in traps , and vor- 
tices in confined geometries 0, . The shell ordering 
in these systems can be altered by the shape or type of 
trap as well as by competing attractive and repulsive in- 
teractions between the particles [H, \^ . In our system, 
the particles within each shell can exhibit an additional 
ordering due to interactions with particles in the neigh- 
boring traps. We show that as a function of substrate 
strength, a rich variety of colloidal cluster crystals can 
be created depending on the filling, and that certain shell 
structures are more stable than others. The number of 
particles in the shells and the number of shells depend on 
the substrate strength, and in certain cases the particles 
form ring structures instead of shells. We also find reen- 
trant ordered phases as a function of substrate strength. 
For weak and strong substrates, the system is ordered; 
however, for substrates of intermediate strength, the sys- 



tem is generally disordered. We find that the depinning 
threshold can show distinct changes at the boundaries 
between these different phases as a function of substrate 
strength, and remarkably, we find that in some cases it is 
possible for the depinning force to decrease with increas- 
ing substrate strength. The velocity-force curves contain 
clear signatures of the different phases corresponding to 
different modes of depinning. 

II. SIMULATION 

We consider a 2D system with periodic boundary con- 
ditions in the x and y directions. The sample contains 
Nc particles interacting with a triangular sinusoidal sub- 
strate containing Ns potential minima. We focus on the 
case where the number of colloids Nc is an integer multi- 
ple of Ns such that / = N^/Ng = n, where / is the filling 
factor and n is an integer. A single particle i responds to 
forces according to the following overdamped equation of 
motion: 

— =Ff + F^ + F'=^*-|-rf. (1) 

The particle-particle interaction force F^"^ — 
— VF(i?ij) and the particle interactions are 
of a Yukawa form, V{Rij) — EQexp{—KRij)/Rij, where 
Eq ~ Z'^* /Aneeoao, e is the solvent dielectric constant, 
Z* is the effective charge, 1/k is the screening length, 
flo is the unit of length which is of order a micron, and 
Rij — ]Ri — Rj]. We model the substrate as a triangular 
lattice, 

F^ = V A sin f ^) [cos(0fe)x - sm{Ok)y] (2) 

where bk = xcos{9k) — ysui(9k) + ao/2, 9i — 7r/6, 
6*2 = 7r/2, and 6^ = 5tt/6. A is the amplitude of the 
substrate potential. The initial colloid positions are ob- 
tained by simulated annealing, where the colloids start 
in a high temperature molten state and are gradually 
cooled to a finite temperature or to T = 0. We have 
also considered other initialization procedures, includ- 
ing placing a commensurate number of colloids in each 
substrate minima and letting the system relax. This 
procedure gives very similar initial states as the simu- 
lated annealing, especially for strong substrates; however, 
for weaker substrates the configuration obtained through 
simulated annealing is generally more disordered, as we 
discuss in more detail later. After annealing, we inves- 
tigate the transport properties by applying an external 
drive F*^^* = Fjjk. We increment Fjj from Fu = to 
some final value in steps of SFjj, and we wait 10^ simula- 
tion time steps between increments to avoid transient ef- 
fects. The depinning threshold is obtained by measuring 
the average colloid velocity V = N''^ J2i {dKi/dt) ■ x. 
The value of SFjj is modified according to the strength 
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of the substrate. For weaker substrates we use smaller 
values of SFjj in order to obtain an accurate depinning 
threshold. 



III. CONFIGURATIONS AND DYNAMICS FOR 

/ = 3, 4, AND 5 

We first consider the cases / = 3, 4, and 5. Earlier 
numerical studies addressed the statics and dynamics for 
/ = 2 and 3 and showed that for strong triangular sub- 
strates, f — 2 produced a herringbone state that transi- 
tioned under an applied drive into an aligned ferromag- 
netic state where the particles move in ID channels ,27J . 
For weaker substrates a triangular lattice formed and the 
depinning was elastic, while for intermediate substrate 
strengths the depinning was plastic and the herringbone 
state depinned into a fluctuating state that reordered at 
higher drives. There was a strong increase in the crit- 
ical depinning force between the elastic and plastic de- 
pinning regimes similar to the peak effect phenomenon 
found at the transition between elastic and plastic depin- 
ning for vortices in type-II superconductors. For / = 3 
the pinned state was always triangular and for weak sub- 
strates the particles depinned into either a moving crys- 
tal or a moving smectic state. At these low fillings it was 
also shown that for weak substrates the critical depinning 
force Fc oa A^, as expected for elastic depinning, while for 
strong substrates, Fc oc A, as expected for single particle 
depinning. 

We now consider several features at / = 3 that were 
not reported in the earlier studies, including the velocity- 
force curves and an additional dynamical phase. Fig- 
ures [Tfa) and [Tfb) show the particle configurations at 
/ = 3 for A — 1.0 and A — 3.0, respectively, while 
Figs.[TJc) and[T{d) show the corresponding structure fac- 
tor S{k). For A = 1.0 in Fig. [TJa,c), a triangular lattice 
forms with one dominant length scale that is the dis- 
tance between the particles. For A = 3.0 in Fig. [2b, d), 
each substrate minimum clearly captures three particles 
and the structure factor has the signature of two length 
scales. There is a well defined hexagonal pattern of max- 
ima in S{k) surrounding the origin in Fig. [Ijd) that was 
not present in Fig. [TJc) which is associated with the 
longer length scale of the substrate minima locations. 
The hexagonal void arrangement at larger k values in 
Fig. [^d) is associated with the smaller length scale of the 
trimer colloid arrangement within each substrate mini- 
mum. In general, as A increases further, the structure 
of the lattice remains the same as shown in Fig. Wic,d) 
except that the trimers gradually reduce in size. 

In Fig. [2] we plot the average particle velocity V ver- 
sus the external drive Fd for A — 1, 2, 3, 4, 5, and 6. 
For the A = 4.0 curve we have labeled the different dy- 
namical phases. The pinned phase (P) transitions into 
a rapidly fluctuating phase we term the random phase 
(R). The random phase is characterized by rapid mixing 
of the particles, as illustrated in Fig. ^a) for A = 4.0 
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FIG. 1: (a, b) The particle locations (black dots) for a tri- 
angular substrate at a filling of / = 3 for different substrate 
strengths of (a) A — 1.0 and (b) A — 3.0. In (b) the trimer 
states that form at each trap are aligned, (c) The structure 
factor S{k) corresponding to panel (a) at A = 1.0. (d) S{k) 
corresponding to panel (b) a,t A — 3.0 shows long range or- 
der and some features at larger k values corresponding to the 
smaller real space length scale of the trimers. 



at Fd — 0.9, where the particle trajectories show signif- 
icant excursions in the direction transverse to the drive 
and are continuously changing with time. In Fig. Eljb), 
the corresponding S{k) taken from a snapshot of the par- 
ticle positions in the random phase shows features of an 
anisotropic liquid. There is triangular ordering in S{k) 
at smaller fc values due to the underlying triangular sub- 
strate; however, the higher order peaks of Fig. [IJd) that 
indicated long range order in the pinned state are re- 
placed in Fig. [3jb) by a ring structure characteristic of 
a liquid. We also find a larger number of peaks along ky 
compared to kx due to the anisotropy induced by the a:- 
direction drive. For higher drives the system transitions 
into a state with no transverse diffusion, as illustrated 
in Fig. [3l[b,d) for a sample with A — 4.0 at Fjj = 1.5. 
In this state, which we term the moving smectic (MS) 
phase, the motion is confined to winding quasi-lD fiows 
as shown in Fig. M.^). The plot of S{k) in Fig. M.d) in- 
dicates that the MS state has stronger partial ordering 
than the random phase but that there is still a tendency 
for the system to be more ordered along the ky direction. 

The R to MS transition is correlated with a decrease 
in V with increasing Fjj, producing a negative dV/dFo 
signature known as negative differential conductivity 
(NDC). Simulations and experiments on superconducting 
vortices moving over periodic pinning sites have revealed 
NDC at transitions from random or turbulent fiows to 
ordered ID flows [i^. The NDC in the vortex system 
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FIG. 2: The average velocity V vs external drive Fd for the 
/ = 3 system in Fig. [l]for A = 1, 2, 3, 4, 5, and 6, from left to 
right. For A = 4.0 the three different phases are labeled for 
the pinned (P), random flow (R), and moving smectic (MS) 
phases. For the A = 2.0 curve there is no random phase but 
two different elastic flow phases EFl and EF2 occur which 
are labeled with lines. 
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FIG. 3: (a) Particle positions (dots) and trajectories (lines) 
over a fixed time interval for the random phase at Fd =0.9 
and A = 4.0 for the / = 3 sample from Fig. [2] The tra- 
jectories are rapidly changing with time and mixing of the 
particles occurs, (c) The corresponding S{k) shows sixfold 
peaks at small k due to the overall triangular ordering of the 
particles on the substrate as well as a ring structure at larger 
k due to the liquidlike nature of the particle positions, (b) 
Particle positions (dots) and trajectories (lines) over a fixed 
time interval in the moving smectic phase at Fd = 1.5 for the 
same sample, where the particles follow well-defined winding 
quasi-lD paths, (d) The corresponding S{k) shows more or- 
dering, including more pronounced peaks along ky that are 
consistent with a smectic phase. 
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FIG. 4: Particle positions (dots) and trajectories (lines) over 
a fixed time interval in the elastic flow phases from the / = 3 
sample in Fig. [2] at A = 2.0. (a) In phase EFl at Fd = 0.25, 
the particles form a triangular lattice and flow through wind- 
ing channels but there is no particle mixing, (b) In phase 
EF2 at Fd = 0.5, the trajectories become straight. The tran- 
sition between these two phases appears as a kink in V near 
Fd = 0.32 in Fig. H 



is much more pronounced than the NDC we observe in 
Fig. [21 and this is likely due to the difference in the type 
of substrate used in the two systems. For the vortex sys- 
tem, artificially fabricated pinning sites produce a short 
range mufhn-tin potential that permits some vortices to 
sit completely outside of the pinning sites in the flat in- 
terstitial regions. In the random phase, a large number 
of vortices move through the interstitial areas and do not 
interact directly with the pins, but at the transition to 
ID flow, these vortices suddenly fall into flowing channels 
that pass through the pinning sites, abruptly increasing 
the effective amount of drag exerted by the pinning sites 
on the vortices. In the sinusoidal triangular substrate we 
consider here, there are no interstitial regions and the 
moving particles are always experiencing some drag from 
the pinning. Additionally, in the vortex system NDC was 
only associated with fillings just above / = 1, the first 
matching field [48], while in Fig. [2] we find no NDC until 
/ > 3. We note that the NDC for the triangular substrate 
is much more robust than NDC in the vortex system since 
it appears for a much wider range of fillings, including 
incommensurate fillings. In a recent vortex experiment 
[49i] NDC did not occur until / > 2.5; however, this effect 
could not be reproduced in vortex simulations using only 
muffin-tin type pinning potentials. This could mean that 
for the particular system considered in Ref. [i^, despite 
the fact that the pinning sites are localized, longer range 
interactions between the pins or distortions of the sub- 
strate may have caused the substrate to behave more like 
the sinusoidal type of substrate we consider here rather 
than like a muffin-tin potential. 

For A > 4.0 we find the same three phases described 
for A = 4.0, with the transitions between phases shifting 
to higher Fu with increasing A. In Fig. [2] we identify 
two distinct phases observed for a weaker substrate with 
A = 2.0. Here the random phase is absent and the system 
depins elastically, with no sharp jump in V at depinning 
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as found for A > 2.0 when the flow above depinning is 
plastic. Even though the depinning transition for ^ = 2.0 
is elastic, we still find transitions between distinct dy- 
namical phases that can be identified through signatures 
in the velocity-force curve. In Fig. [2] we mark the two 
elastic flow phases EFI and EF2 that appear on cither 
side of a discontinuity in dV/dFu centered at Fo = 0.32. 
In both phases the particle lattice is triangular and has 
an S{k) similar to that shown in Fig.[ljb). The EFI and 
EF2 phases also appear for A — 3.0 with the kink in V 
shifted up to higher . In Fig. Ufa) we plot the particle 
trajectories in phase EFI for the A — 2.0 system. The 
triangular particle lattice moves in a zig-zag pattern in 
order to avoid passing over the potential maxima of the 
triangular substrate. At higher drives the particles no 
longer avoid the potential maxima and enter phase EF2 
where they move in the direction of drive, as illustrated 
in Fig.SJb). 

For / = 4, when the substrate is weak a triangular 
lattice of particles forms, as illustrated in Fig. [5Ua,b) for 
a sample with A = l.Q. Weak substrates still affect the 
lattice by breaking rotational symmetry and causing the 
particle lattice to preferentially orient in a direction de- 
termined by the filling factor. At / = 4, one of the 
symmetry axes of the particle lattice is aligned with the 
X axis as in Fig. [5]^ a), while in Fig. [TJa) at / = 3 it 
was aligned with the y axis. As the substrate strength 
increases for / = 4, quadrimer arrangements of the parti- 
cles form in each substrate minimum, and the orientation 
of neighboring quadrimers varies with A. For A = 8.0, 
Fig. [SJc) shows that all the quadrimers are aligned in 
the same direction. The corresponding S{k) in Fig. [Sl^d) 
shows sixfold peaks at small k from the triangular sub- 
strate, while at larger k a fourfold void structure ap- 
pears due to the square ordering of each quadrimer. At 
A = 11 in Fig. [DJe), the individual quadrimers remain 
intact but their long range ferromagnetic orientational 
ordering is lost, producing a ring structure in S(k) as 
shown in Fig. [5ff). We have tried several different meth- 
ods of preparing the initial configuration for / = 4 and 
A = 11.0 but have not found a ferromagnetically ordered 
state. 

The orientational ordering of dimer and trimer states 
was previously shown to arise due to quadrupole or higher 
pole moment interactions between neighboring rt-mers. 
For quadrimers, it is likely that the higher pole moments 
play a more important role. If the pole moment favors 
ferromagnetic alignment, than as A increases the size of 
the quadrimer decreases, reducing the pole moment re- 
sponsible for the orientational ordering. Previous stud- 
ies foimd that as the substrate strength increases, the 
temperature at which dimers and trimers lose their ori- 
entational ordering drops due to the decreasing effective 
pole moment. We note that the configurations in Fig. [5] 
are all obtained at T = 0.0. When we anneal the system 
from a high temperature and slowly decrease the temper- 
ature to zero, we still obtain the same states illustrated 
in Fig. [HI suggesting that the energy of the ferromagnetic 
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FIG. 5: (a,c,e) The particle configurations and (b,d,f) cor- 
responding S{k) for samples with / = 4. (a,b) At A = 1.0 
the particles form a triangular lattice. (c,d) At A = 8.0 a 
quadrimer state forms. The quadrimers have a ferromagnetic 
alignment. (e,f) At A = 11.0 the quadrimer state has no 
orientational ordering, producing ring structures in S{k). 



state must be only slightly smaller than that of the ro- 
tationally disordered state. As A increases, the picture 
of rigid quadrimers begins to break down, as shown in 
Fig.EJe) where the quadrimers become slightly distorted. 
This could further decrease the contributions of differ- 
ent pole moments to the orientational ordering or may 
even lead to competing interactions between neighbor- 
ing quadrimers, resulting in a ground state with ordering 
produced by very long range interactions, as found in cer- 
tain spin ice systems. For A > 11.0 we find states similar 
to those shown in Fig.[5je,f) with increasing distortion in 
the square ordering, until for high enough A the distor- 
tion becomes strong enough to permit some particles to 
sit at the center of a substrate minimum. The dynamic 
phases for / = 4 are similar to those illustrated in Fig. [5] 
for the / = 3 system, with elastic depinning for low A 
crossing over to plastic depinning at higher A, followed 
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by a dynamical reordering transition at higher Fd. 

For / = 5 at Fd = 0, as a function of A we find disor- 
dered phases interspersed among ordered phases. This is 
also the first filling at which the n-mer assumption breaks 
down and a shell structure of the particles in each sub- 
strate minimum begins to form, in contrast to the case 
of isolated trapped clusters of repulsively interacting par- 
ticles, which first develop equilibrium shell structure at 
f — 6 [35l.[37ll43l. [50l - l53j but have a metastable state with 
one particle in the center for / = 5 [HJ] . In Fig. [D^a) for 
yl = 0.5at/ = 5,a triangular lattice containing a small 
tilt distortion forms. The smearing of the correspond- 
ing S{k) in Fig. MJd) indicates that for weak substrates 
at this filling, the particle lattice is not commensurate 
with the substrate. At A = 1.75, shown in Fig. [ljc,d), 
the lattice is disordered as indicated by the ring pattern 
in S{k). When we vary the initialization protocol, we 
always observe a disordered lattice for 0.9 < A < 2.0. 
For 2.0 < A < 7.0 the system orders again as shown in 
Fig.[6lje,f) for A = 5.0. Each substrate minimum contains 
what we term a jack pattern consisting of one particle lo- 
cated at the center of the minimum surrounded by four 
particles in a square arrangement. As indicated by the 
dashed lines in Fig. Elje), the jacks are tilted at -|-20° or 
—20° to the X axis in every other row, producing a her- 
ringbone structure. The ordering is not complete as there 
are local regions where the jack structure breaks down, 
producing some smearing in S{k) as plotted in Fig. [6][f). 
The structure can be viewed as the beginning of a shell 
ordering where, in each substrate minimum, the first shell 
contains a single particle and the next shell contains four 
particles. For A > 7.0 the interactions between parti- 
cles in neighboring substrate minima are reduced and the 
pentagon ground state structure expected for particles in 
an isolated trap emerges. The pentagons exhibit ferro- 
magnetic order and are all aligned in the same direction, 
as illustrated in Fig. |njg,h) for A — 11.0. For increasing 
A the pentagons become smaller; however, we observe no 
further change in the structure. 

We can relate the different pinned structures at / = 5 
to features in the velocity force curves and to dynamic 
phases. For A < 1.5 there is an elastic depinning into a 
moving triangular lattice. For 1.5 < A < 2.4 the disor- 
dered phase depins into a fluctuating random phase sim- 
ilar to that observed for / = 3, and at higher drives the 
system can reorder into a moving crystal (MC) phase. 
In Fig. [Tlja) we plot V versus Fjj for a sample with 
A — 2.15. The random flow phase is characterized by 
large fiuctuations, and is followed by a sharp transition 
near Fjj = 0.034 into a MC phase with small fiuctua- 
tions. The transition between these two phases can also 
be detected by analyzing the noise fiuctuations of the 
velocity response. In the random phase the fiuctuations 
are characterized by a broad band noise signal with 1//" 
characteristics where a — 1.5 to 2.0, while in the MC 
phase there is a narrow band signal with a character- 
istic frequency that increases with increasing Fjj., simi- 
lar to what has been observed for dynamically reordered 
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FIG. 6: (a,c,e,g) The particle configurations and (b,d,f,h) 
corresponding S{k) for samples with / = 5. (a,b) At A = 0.5 
the particles have triangular ordering that is not commen- 
surate with the substrate, as indicated by the smearing of 
the peaks in S{k). (c,d) At A = 1.75 the system is disor- 
dered. (e,f) At y4 = 4.1 the system forms what we term a jack 
state with one particle located at the center of each potential 
minimum surrounded by four outer particles in a square con- 
figuration. In large regions the squares form a herringbone 
structure, as indicated by the dashed lines; however, some lo- 
calized disorder is present in the sample. (g,h) At A = 11.0, 
an aligned pentagon state with long range orientational or- 
der forms. For A > 11.0 the same structure persists but the 
pentagons shrink in size. 
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FIG. 7: The velocity V vs force Fo curves for samples with 
/ = 5. (a) At A — 2.15 we observe a pinned phase (P), a 
random phase (R), and a moving crystal phase (MC). (b) At 
A = 4.1 the system depins into a moving channel phase (C) 
illustrated in Fig. [8ja). This is followed by the R phase and 
the moving smectic phase (MS), (c) At A = 7.2 the C phase is 
absent and the system transitions directly from the R phase 
to the MS phase. 




(a) 




(b) 



FIG. 8: Particle positions (dots) and trajectories (lines) over 
a fixed time interval for the / = 5 sample, (a) The moving 
channel (C) phase at ^4 = 4.1 and Fd ~ 0.05. (b) The moving 
smectic (MS) phase at A = 7.2 and Fd = 1.1. 



phases of superconducting vortices moving over periodic 
substrates JSl- For 2.5 < A < 3.3 we find no dynamical 
reordering and the particle arrangement remains disor- 
dered up to the highest drives we considered; however, it 
is possible that a transition into a moving ordered state 
could occur at very high drives. For 3.1 < A < 6.1, where 
the pinned system is in the ordered jack state, the depin- 
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FIG. 9: (a) The dynamical phase diagram Fd vs A for a 
sample with / = 5. Regions where the P, C, MC, R, and 
MS phases occur are marked, (b) The depinning curve Fc 
from panel (a) plotted on a log-log scale. Dashed lines mark 
the transitions between the pinned triangular crystal (PC), 
pinned disordered or pinned glass (PC), pinned jack (PJ), 
and pinned pentagon (PP) states, (c) The critical depinning 
force Fc vs A for / = 1, 3, 4, and 5, from top to bottom. 
Upper dashed line: a linear fit Fc oc A. Lower dashed line: 
a fit to Fc oc A^. These curves indicate that, in general, the 
system exhibits collective depinning behavior. 



ning is elastic and the particles depin into the ordered 
winding channel flow phase (C) illustrated in Fig. [Sja). 
As F]j increases, there is a crossover to a random phase 
followed by another transition at high drives into a mov- 
ing smectic phase (MS) of the type shown in Fig. [Hl^b). 
The velocity signatures associated with these transitions 
appear in Fig. [7]^b) for a sample with A — 4.1. The C 
phase is associated with small velocity fluctuations and 
is followed first by the strongly fiuctuating R phase and 
then by the MS phase at higher Fd- For A > 6.0 the 
C phase is lost and the sample depins from the ordered 
pentagon state into a random flow state, as illustrated in 
Fig. \T[^c) for a sample with A = 7.2. This is followed by 
the formation of a moving smectic state at higher Fd. 
As A increases, the extent of both the pinned region and 
the random fluctuating state increases. 

By conducting a series of simulations, we map the dy- 
namic phase diagram shown in Fig.lHJa). The disordered 
R phase separates the moving crystal and moving smectic 
phases. There is also a domelike region where the wind- 
ing channel C phase occurs. In Fig. Mjo) we plot the de- 
pinning line on a log-log scale with dashed lines marking 
the regions where the pinned phases illustrated in Fig. [S] 
occur. For A < 0.9, the depinning force Fc monotonically 
increases with increasing A. There is a small decrease in 
Fc at the onset of the pinned disordered or pinned glass 
(PG) phase, followed by a plateau in Fc. Near A = 3.0, 
Fc begins to increase rapidly with increasing A in the 
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pinned jack state, while in the pinned pentagon state 
Fc stih increases with increasing A but with a reduced 
slope. The behavior of Fc near the onset of the pinned 
disordered phase is somewhat unusual since it indicates 
that even though the substrate strength is increasing, 
the depinning force does not increase. When the system 
is disordered, there are numerous dislocations present in 
the particle structure that produce weak spots which flow 
first at the depinning transition. In contrast, in the de- 
pinning of a crystalline state, there are no weak spots 
since the lattice structure is the same everywhere. This 
behavior is the opposite of the so-called peak effect phe- 
nomenon found for vortices in type-II superconductors, 
where a sudden increase in Fc occurs when the vortex 
lattice becomes disordered. The peak effect is generally 
believed to arise due to the softening of the vortex lat- 
tice when dislocations are present. In the softer lattice, 
the vortices can more freely move to adjust to a ran- 
dom pinning potential and maximize the pinning force. 
In our system, the pinning potential is not random, so 
the crystalline states are more strongly pinned than the 
disordered states. This has been demonstrated clearly in 
studies of commensurate-incommensurate transitions for 
both vortices and colloids in periodic substrates, where 
crystalline states form at integer values of /. Near but 
not at commensuration, numerous defects appear in the 
commensurate lattice, and these defects cause a reduc- 
tion in Fc- As a result, Fc passes through a series of peaks 
at integer values of / as the filling fraction is varied. In 
the system we are considering here, the / = 5 state is 
unusual since even through the system is at an integer 
filling, we find a regime where the pinned state is non- 
crystalline. This is in contrast to the lower integer fillings 
/ = 1, 2, 3, and 4, where all the pinned states are or- 
dered. This disordering effect and the plateau or drop in 
Fc at intermediate substrate strengths in the disordered 
pinned regions is a general feature of / > 5 systems at 
integer fillings, and both features becomes more promi- 
nent for higher fillings. It may appear from Fig.[SJb) that 
the plateau in Fc occurs in the disordered region and not 
in the pinned jack state; however, the plateau persists 
into a window of the pinned jack state since the PJ state 
contains some defects. 

In Fig. ^c) we plot Fc vs A for samples with f — 1, 
3, 4, and 5. The dynamics and configurations of the 
/ = 2 system were described in detail in previous work 
27| . For / = 1 we find that Fc oc A, as indicated by the 
upper dashed line. This is expected in the limit of single 
particle depinning and indicates that particle-particle in- 
teractions are unimportant at depinning for this filling. 
For the higher fillings / > 1, we instead find Fc oc A^ 
as indicated by the lower dashed line. This is character- 
istic of collective depinning transitions, where the inter- 
actions among the particles within each minimum play 
an important role in the depinning process. In the limit 
of very strong substrates, the depinning generally occurs 
when one particle is pushed near the saddle point of the 
substrate by the other particles in the potential mini- 
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FIG. 10: Fc versus A. (a) At / = 6, a sharp increase in 
Fc occurs at the transition from a pinned triangular lattice 
to a disordered state. Dashed lines indicate regions where 
different pinned phases appear. Triangular: triangular lat- 
tice; Disordered: disordered state; Pentagon/Triangular: a 
monomer at the center of each well surrounded by five par- 
ticles arranged in a pentagon; Hexagon: a hexagonal ring 
in each substrate minimum, (b) At / = 7, the triangular 
to disordered transition is associated with a small dip and a 
plateau in Fc with increasing A. Dashed lines indicate the re- 
gions where different pinned states appear. Triangular: trian- 
gular lattice; Disordered: disordered state; Pentagon/Dimer: 
a composite pentagon-dimer arrangement in each minimum; 
Hexagon/Monomer: a composite hexagon- monomer arrange- 
ment in each minimum; Ring: a single ring of particles in each 
minimum found at large A. 



mum. This single particle suddenly escapes and triggers 
a cascade of escape events in the rest of the system. At 
intermediate substrate strengths, the higher order fillings 
all generally show a change in the slope of Fc versus A 
associated with a change in the n-mer structure. 



IV. CONFIGURATIONS AND DYNAMICS FOR 

6 < / < 10 

In Fig. rTOT a) we plot Fc versus A for a sample with 
/ = 6. A large jump in Fc occurs a.t A = 2.5 at the 
transition from a pinned triangular lattice, illustrated in 
Fig. irTTa). to a disordered pinned state, illustrated in 
Fig. fTlTb). For 2.5 < 5.5, the system becomes disor- 
dered, as shown in Fig. ITlT b) for A = 2.75, and has a 
ringlike S{k) signature, shown in Fig. [TTJg). In the dis- 
ordered region the system depins into a plastic random 
flow state, while for A < 2.5 the system depins elastically 
into a MC state. At this fllling, Fc is depressed in the 
triangular state since it is not possible for all the parti- 
cles to simultaneously be in a triangular lattice and sit at 
substrate minima locations, so that some particles are in- 
stead located on substrate maxima. When the substrate 
strength increases enough, the particles are forced off of 
the substrate maxima and the system disorders, allowing 
more particles to be located closer to substrate minima. 
For 5.5 < ^ < 15, a shell structure emerges. The sys- 
tem forms a pentagonal/triangular lattice composite in 
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(a) X (b) X (c) x_ 




FIG. 11: The particle positions after annealing for a sample 
with / = 6. In panels (c,d,e), dashed lines indicate the ori- 
entation of the pentagons in selected traps, (a) At A — 1.0 
and Fd = a triangular state forms, (b) At A — 2.75 and 
Fd = a disordered state occurs, (c) At j4 = 6.0 and Fd ~ 
there is a pentagon/monomer state with no orientational or- 
der, (d) At A — 12.0 and Fd = we find an aligned pen- 
tagon/monomer state, (e) At A — 12.0 and Fd = 0.45, just 
below the depinning transition, the pentagon/monomer state 
is polarized by the applied drive, (f) At A = 17.0 and Fd = 0, 
a hexagonal ring state with orientational order forms. Dashed 
lines indicate the two hexagon alignments in representative 
traps, (g) S{k) for the disordered phase in panel (b). (h) 
S{k) for the hexagonal ring state in panel (f). 



which one particle sits at each substrate minimum and 
is surrounded by five particles in a pentagon arrange- 
ment, as illustrated in Fig. [TIJc) for A = 6.0. In this 
case the ordering is not ferromagnetic and the tilt of the 
pentagons varies from site to site; however, for A > 8.0 
the pentagon-monomer states become aligned, as shown 
in Fig. [TlTd) for A = 10. An apphed drive can induce a 
polarization of the pentagons in the driving direction, as 
illustrated in Fig. \TT\ e] for A = 12.0 under an x direc- 
tion drive just before depinning. For yl > 15 a ring state 
forms in which the center of the substrate minimum is 
empty and the particles form a hexagon. The rings are 
not all aligned; instead, two orientations coexist in the 
sample as shown in Fig. Illf f). This produces twelve- 
fold modulations in the structure factor as illustrated in 
Fig. fTTT h). As A increases further, the ring structure per- 
sists but decreases in diameter. In general we find ring 
structures rather than shell structures for / > 6 at large 
A. The different phases are highlighted in Fig. fTUfa). 
Figure [TU^b) shows that a similar set of pinned phases 
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FIG. 12: The particle positions after annealing for a sample 
with / = 7. (a) The disordered state at A = 4.0. (b)At 
A = 12.0 a pentagon/dimer composite appears in which 
the dimers are partially aligned. Dashed lines indicate two 
different pentagon/dimer orientations, (c) At A = 20.0 a 
hexagon/monomer composite forms. Dashed lines indicate 
two different hexagon orientations, (d) At A — 25.0 a ring 
structure appears with seven particles per ring. 



forms in a sample with f — 7. There is still a triangular 
lattice for small A; however, sd A = 3.75 when the system 
enters the disordered phase, Fc shows a small decrease 
followed by a plateau up to the end of the disordered 
phase at ^ = 7.0. For / = 7 the triangular lattice is 
more commensurate with the substrate and thus better 
pinned than the triangular lattice at low A in the / = 6 
system. In the disordered phase, illustrated in Fig. [T2la) 
for A = 4.0, S{k) contains a ringlike structure similar 
to that in Fig. \lV[g). For 7.0 < A < 16 the system 
forms a shell structure with an outer pentagon and two 
dimerized inner particles, as shown in Fig. I12r b) for A — 
12.0. The pentagon/dimer structures are only locally 
aligned, leading to smearing in the corresponding S{k). 
Figure [I2\ c] illustrates the configuration at A = 20.0 
where an array of hexagons appears with a monomer at 
the center of each hexagon. The hexagons have one of two 
possible orientations, similar to the / = 6 hexagon state 
shown in Fig. [THf). For A = 25 shown in Fig. \W[d), 
the shell ordering is replaced by a ring structure where 
each ring has seven particles and the ring radius decreases 
with increasing A. As also found for the / = 6 case, it 
is possible for an external drive to orient the structures 
close to depinning. 

The dynamical phases for / = 6 and f — 7 are very 
similar to those observed at / = 5. The system depins 
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FIG. 13: The particle positions after annealing for a sample 
with / = 8. In (a,b), dashed circles indicate the centers of 
the substrate minima, (a) At ^4 = 2.0 a triangular lattice 
forms that is rotated by 19.1° with respect to the underlying 
substrate. The dashed line indicates the direction along which 
8 particles fit between consecutive substrate minima, (b) At 
A = 5.0, dimerization of a subset of the particles occurs and 
the dimers are aligned approximately —41° from the x axis. 
The dashed line indicates the lattice alignment direction, and 
a representative hexagon surrounding one dimer is shown, (c) 
At A = 20.0, seven particles surround a single monomer in 
each minimum, (d) At j4 = 26.0 a ring state appears with 
eight particles per ring. 



elastically at small A in the triangular ordering regime 
and depins plastically in the disordered pinned regime. 
A disordered flow regime extends to divergingly large 
drives for regions of A where a disordered pinned phase 
is present. For values of A above the disordered pinned 
phase, the strongly driven system dynamically orders into 
a moving smectic phase, and the driving force at which 
this transition occurs increases with increasing A. For 
the higher fillings of / = 6 and / = 7, there are some 
additional features in the V — Fn curves at high A that 
did not appear at / = 5. For example, in the fluctuating 
phase, the system shows a coexistence regime with some 
particles that remain pinned while other particles move 
along ID channels or troughs without rearranging the 
pinned particles. This feature becomes more noticeable 
at higher /. 

For f — S there are several new features in the dynamic 
phases due to the fact that there is no intermediate dis- 
ordered pinned phase. This results from a matching of 
the registry of the f — 8 lattice with the underlying sub- 
strate. At A = 2.0 in Fig. [TSf al. the particles form a 
triangular lattice that is slightly anisotropic and rotated 




A 



FIG. 14: The dynamical phase diagram Fd vs A for a sam- 
ple with f — 8. The pinned (P) phases include: a triangu- 
lar lattice (Triangular), a triangular lattice containing dimers 
(Triangular/Dimer), heptagons surrounding monomers (Hep- 
tagon/Monomer), and a ring state (Ring). There is a large 
increase in Fc at the onset of the dimerized state. The tran- 
sition flow region is labeled Tr. The MC, MS, and R flow 
phases are also marked. 



by 19.1° with respect to the underlying substrate. This 
orientation permits as many particles as possible to sit 
close to substrate minima, and as indicated by the dashed 
line in Fig. I13f a). along one particle lattice symmetry di- 
rection, there are eight particles filling the space between 
consecutive substrate minima. This arrangement differs 
from the disordered lattice found at the eighth match- 
ing field in a sample with a muffin-tin potential [55J . As 
A increases, each substrate minimum traps two particles 
that dimerize and are surrounded by six other particles, 
as illustrated in Fig. [TSf b). The dimers are rotation- 
ally ordered, and the overall particle lattice is rotated by 
40.89° with respect to the substrate, as indicated by the 
dashed line in Fig. [T^ b) . Orientationally ordered dimers 
did not form for / = 5, 6, or 7. For A > 17.0 the dimer 
state is lost and replaced by a monomer at each substrate 
minimum surrounded by an outer shell of seven particles, 
as shown in Fig. ^c) for A = 20.0, and for A > 25.0 
a ring state with eight particles per ring occurs, as illus- 
trated in Fig.[n);d) for A = 26.0. 

The lack of disorder for / = 8 is clearly shown in the 
dynamical phase diagram in Fig. 1141 There are no dis- 
ordered pinned phases, and the random flow phase (R) 
does not extend out to large Fd but is instead bounded 
on the high Fd side by the MS flow state. For A < 2.0, 
the triangular lattice depins elastically in the direction of 
drive into an MC state. There is a sharp increase in Fc at 
the onset of the dimerized state, and the MC flow is re- 
placed by MS flow at higher drives. For A < A < 6.0, the 
system depins transverse to the driving direction along 
one of the symmetry directions of the particle lattice at 
—40.89° from the x axis, as illustrated in the inset of 
Fig. [ini As Fd increases, there is a sharp transition out 
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FIG. 15: Main panel: Vx (upper curve) and Vy (lower curve) 
vs Fd for a sample with / = 8atA = 5.0. P: pinned phase; 
Tr: transverse flow phase; R: random flow phase with motion 
oriented along the x direction only. There is a sharp transition 
between the Tr and R flow phases. Inset: Particle positions 
(dots) and trajectories (lines) over a flxed time interval for the 
same sample in the transverse flow phase just above depinning 
at Fd ~ 0.0125. Dashed circles indicate the centers of the 
substrate minima. The particles do not flow with the applied 
drive along the a;-direction but instead travel at —40.89° to 
the X axis along one of the principal axes of the particle lattice. 



of the transverse flow regime into a random flow phase 
with an average velocity oriented in the driving direc- 
tion, as shown in the main panel of Fig. 1151 where we plot 
(Vx) and (Vy) versus Fd. In the transverse flow regime 
for 0.012 < Fd < 0.0175, {Vy) is negative and (T4) is 
positive, and both and (Vy) increase Unearly with 

increasing drive until, at Fd = 0.0175, (Vy) drops to 
zero and (Vx) jumps up at the transition to the random 
flow regime. Note that the orientation of the pinned lat- 
tice is degenerate, so for samples prepared with different 
random seeds, the particles would flow along —40.89° in 
half the cases and along -1-40.89° in the other half of the 
cases. The transverse flow regime forms an intermedi- 
ate state between the pinned and random flow regimes. 
Flows where the particles do not move in the direction of 
drive but at an angle to the drive have previously been re- 
ported for colloids on triangular substrates at a filling of 
/ = 2.0 [5^. In that study, the drive was rotated 90° rel- 
ative to the substrate symmetry direction from the case 
considered here. Dimerization of all the particles in the 
system created an orientational degree of freedom that 
determined the flow direction of the colloids just above 
depinning. Numerical studies of vortices in type-II super- 
conductors with kagome and honeycomb pinning arrays 
also produced similar flows at angles to the driving di- 
rection in cases where two vortices were located in the 
center of the kagome or honeycomb plaquettes, forming 
dimer states 57 [. 

9, triangular configurations ap- 
2.0. This is 




In samples with / 
pear at low A, as shown in Fig.fTCTa) for A 



FIG. 16: (a-f) The particle positions after annealing, (a) 
For / = 9 at ^ = 2.0, a triangular lattice forms, (b) For 
/ = 9 at = 7.0, each substrate minimum contains a linear 
trimer state aligned along the x direction surrounded by dimer 
states, as indicated by the dashed lines, (c) For / = 9 at 
A — 15.0, triangular trimers sit at each substrate minimum 
surrounded by a hexagon of six particles, (d) For / — 10 
at yl = 1.0, a partially ordered triangular state occurs, (e) 
For / = 10 at A = 12.0, aligned quadrimers form at each 
substrate minimum surrounded by six particles, (f) For / — 
10 at j4 = 15.0, aligned triangular trimers appear that are 
surrounded by seven particles, (g) Fc vs A for / = 9 (upper 
curve) and / = 10 (lower curve). 



the same configuration found at the ninth matching field 
for vortices in triangular pinning arrays in Ref . [ssf . As A 
increases, the particles align with the x axis in an unusual 
new structure illustrated in Fig. [THT b) for A — 7.0. The 
system breaks into dimer and linear trimer states. Each 
substrate minimum contains a linear trimer flanked by 
two dimers, all oriented in the x direction, while between 
adjacent minima there are elongated dimers oriented in 
the y direction. For A = 15.0, Fig. [Wl c) shows that a 
superlattice of aligned triangular trimers appears at each 
substrate minimum, surrounded by an outer shell of six 
particles forming a hexagon. As A increases further, we 
find a transition to a state with a monomer at each sub- 
strate minimum surrounded by an eight-particle shell, 
followed by a ring state at high A (not shown). For the 
parameters we considered, we never observed a state with 
two dimerized particles at the centers of the substrate 
minima surrounded by seven outer shell particles. Either 
this state appears only for an extremely narrow range of 
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A, or it is simply too high in energy to form at alL In 
general we find that at the higher filhngs, certain shell 
combinations are not observed. These are often, but not 
always, associated with incommensurate structures such 
as odd-even shells like the dimer-heptamer structure for 
/ = 9. 

For / = 10 at A = 1.0, Fig. indicates that a 

partially ordered triangular lattice appears. This is a 
floating disordered solid that is only very weakly pinned 
by the substrate. At A = 1.5, the system transitions 
to a pinned disordered solid or pinned glass state (not 
shown) , where the last vestiges of triangular ordering are 
lost. For A « 7, the particles begin to develop an in- 
cipient quadrimer structure around each substrate min- 
imum. The quadrimers are still very extended and the 
lattice is only partially localized by the substrate minima. 
At A = 9.0, the quadrimers become well localized within 
the minima and the state illustrated in Fig. [lE\ e) for 
A = 12.0 emerges, with aligned quadrimers surrounded 
by hexagons. At A — 15.0, aligned triangular trimers sit 
at each substrate minimum as illustrated in Fig. I16f f). 
This is similar to the state shown in Fig. llGf c) except 
each trimer is surrounded by seven particles in the outer 
ring instead of six. A similar set of orderings as those 
found for f — 9 occurs for / = 10 as A increases until 
a ring state forms at high A. In Fig. [TBT g) we plot Fc 
versus A for / = 9 and / = 10. For / = 9, the triangu- 
lar lattice in Fig.fTCTa) transitions into the aligned linear 
trimer state in Fig. [TCTb) at A — 7.0, and there is a kink 
in Fc associated with this transition. For / — 10, the low- 
A, partially ordered lattice illustrated in Fig. [THl d) has a 
very low Fc since it is still floating above the substrate. 
A rapid increase in Fc at A = 1.5 is correlated with the 
transition to a pinned disordered solid. Fc gradually in- 
creases as the quadrimer structure gradually organizes, 
until at A = 9, when the quadrimers become well local- 
ized within the substrate maxima into the state shown in 
Fig. [TST e). Fc kinks upward again. Near A — 20 there is 
another upward kink in Fc at the formation of the ring 
state, which depin as if they have no orientational degree 
of freedom. 



V. HIGHER ORDER FILLINGS AND 
DISCUSSION 

For higher order fillings / > 10, we generally observe 
the same features described for the / > 5 fillings, so we 
believe that we have identified the generic features that 
arise in this type of system. For higher fillings it is likely 
that instead of only two shells of particles, three or more 
particle shell structures will form. The orientational or- 
dering of particles trapped in adjacent substrate minima 
is also likely to become more fragile as the number of 
shells increases and the outer shell becomes more and 
more circular. Thermal effects on multiple shell struc- 
tures are beyond the scope of this work; however, we 
expect that the shell structures would have multiple dis- 
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FIG. 17: The particle configurations after annealing for high 
filling fractions, showing orientationally ordered states, (a) 
For / = 11 at ^ = 14.0, each substrate minimum contains an 
inner shell of five particles and an outer shell of six particles, 
(b) For f = 12 ai A — 1.0, a triangular lattice appears, (c) 
For / = 12 at A = 10.0, there are four particles in the inner 
shells surrounded by eight particles, (d) For / = 12 at j4 = 20, 
there are five particles in the inner shells surrounded by seven 
particles. 



ordering transitions as a function of temperature that 
are distinct from the transitions observed at lower fill- 
ings in colloidal molecular crystals [H, [2^ . Another fac- 
tor that could infiuence the ordering is the shape of the 
substrate minima. Here we considered a sinusoidal sub- 
strate; however, with optical trapping, other forms of 
substrates could be created which may favor or disfavor 
certain shell structures from forming. 

In Fig. [T7] we show some representative features of 
higher filling states. At / = 11 and A = 14.0 in 
Fig. ll7f a). an aligned pentagonal inner ring forms at each 
substrate minimum with an outer hexagonal ring. For 
f — 12, a. triangular lattice appears for A — 1.0 as illus- 
trated in Fig. I17f b). an oriented structure with four par- 
ticles in the inner shell and nine in the outer shell forms 
at A = 10.0 as shown in Fig.IUl^c), and at A ^ 20.0 there 
are five particles in the inner shell and seven in the outer 
shell as seen in Fig. [TTT dV 



VI. SUMMARY 

We have examined the static configurations and driven 
dynamics of particles interacting with a repulsive Yukawa 
potential, such as colloids, in the presence of a triangular 
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substrate for varied integer fillings of up to twelve parti- 
cles per substrate minimum. Under static conditions, we 
observe a rich variety of crystalline structures for fillings 
that have not been previously reported for this model. 
We show that for fillings / > 4, shell structures can form 
in the substrate minima with well-defined numbers of 
particles in the inner and outer shells. These shell clus- 
ters can exhibit orientational ordering across the sample. 
As we vary the substrate strength for a fixed filling level, 
we find that a series of different particle lattices occur 
with diflFerent numbers of particles in the shells and dif- 
ferent amounts of orientational ordering. In the limit of 
strong substrates we observe a ring lattice. Several fill- 
ings exhibit disordered particle configurations; however, 
these fillings can show a reentrant ordering behavior as 
the substrate strength is increased. Under an applied 
drive, the different particle orderings alter the critical 
depinning force, producing either a strong increase in 
the depirming force or a plateau and even decrease in 



the depinning force at structural transitions induced by 
increasing the substrate strength. We also find a remark- 
ably rich variety of dynamical phases that can be corre- 
lated with the structure of the static pinned configura- 
tion. As the drive increases, we find moving crystalline 
statc;s and moving smectic states as well as disordered 
flow states. Our results are relevant to colloids on peri- 
odic substrates, vortices in nanostructured superconduc- 
tors, and other commensurate-incommensurate systems, 
including aspects of friction and nonequilibrium trans- 
port. 
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